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ABSTRACT 

We determine the critical size Zcrit of a filament of cold [T ~ 10'' K) gas that is 
in radiative equilibrium with X-ray emitting gas at temperatures Tout ^ 10^ — 
10^ K. Filaments smaller than Zcrit will be rapidly evaporated, while longer ones 
will induce the condensation of the ambient medium. At fixed pressure P, Icmt 
increases as T^^/^, while at fixed Tout it scales as P^^ . ft scales as /^/^, where / is 
the factor by which the magnetic field depresses the thermal conductivity below 
Spitzer's benchmark value. For plausible values of /, Zcrit is similar to the lengths 
of observed filaments. In a cluster such as Perseus, the value of Zcrit increases by 
over an order of magnitude between the centre and a radius of 100 kpc. If the 
spectrum of seed filament lengths / is strongly falling with /, as is natural, then 
these results explain why filaments are only seen within a few kiloparsecs of the 
centres of clusters, and are not seen in clusters that have no cooling flow. We 
calculate the differential emission measure as a function of temperature for the 
interface between filaments and ambient gas of various temperatures. We discuss 
the implications of our results for the origin of the galaxy luminosity function. 

Key words: conduction - cooling flows - galaxies: clusters: general - galaxies: 
formation 



1 INTRODUCTION 

Intergalactic space within clusters of galaxies contains 
enormous quantities of gas that is at the virial temper- 
ature. In about three quarters of all rich clusters, the 
temperature falls by a factor of 2 to 3 within a few tens 
of kiloparsecs of the cluster centre - such clusters are 
said to possess a 'cooling flow'. Nearby examples are the 
Virgo, Hydra and Perseus clusters. 

Optical emission-line fllaments are observed around 
NGC 1275 at the centre of the Perseus cluster (Minkowski 
1957; Lynds 1970; Conselice, Gallagher & Wyse 2001, 
and reference therein). Similar structures of extended 
Ha emission are seen near the centres of other cooling- 
flow clusters (e.g., Virgo, Hydra, A2597, A2052, A1795; 
see Heckman 1981; Hu, Cowie & Wang 1985; Johnstone, 
Fabian & Nulsen 1987; Heckman et al. 1989). The ori- 
gin of these fllaments and the mechanism responsible for 
their ionization have been extensively studied, but defini- 
tive conclusions have not been reached. Two formation 
scenarios have been proposed. In the first, the fllaments 
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form as the intracluster medium (ICM) cools (e.g., Fabian 
& Nulsen 1977, Cowie, Fabian & Nulsen 1980). In the 
other scenario, fllaments are produced through cosmic 
infall of either cold gas or gas- rich galaxies (e.g., Soker, 
Bregman & Sarazin 1991; Baum 1992; Sparks 1992). The 
fact that filaments are observed only in clusters with 
cooling flows has often been considered evidence for the 
first scenario, but by no means all cooling-fiow clusters 
have detectable optical emission lines (e.g., A2029; see 
Fabian 1994). Some observational data favour the infall 
scenario. For example. Sparks, Macchetto & Golombek 
(1989) showed that filaments contain dust with normal 
Galactic extinction properties, which fllaments formed 
by condensation of the ICM should not (see also Don- 
ahue & Voit 1993). In addition, the fllaments appear to 
be dynamically disturbed systems, suggesting that they 
form through violent processes, such as merging, rather 
through the quiescent condensation of the X-ray emitting 
ICM. 

An important related question is, what mechanism 
is responsible for heating the filaments? Many possible 
ionization sources have been considered: Active Galac- 
tic Nuclei (AGN); shocks; massive stars; radiation from 
the ICM; magnetic reconnection (e.g., Heckman et al. 
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1989). None of these ionization mechanisms, individu- 
ally, can satisfactory account for the observed emission- 
line ratios and intensities, so possible combinations of 
these processes have been explored (e.g., Sabra, Shields 
& Filippenko 2000). However, there are indications that 
the heating source is physically associated with the fil- 
aments (Conselice et. al 2001 and references therein). 
Thus, it is likely that stellar photoionization makes an im- 
portant contribution to powering their optical emission- 
lines, since young stars are systematically found in the 
filaments (e.g., McNamara, O'Connell & Sarazin 1996). 

The existence of cool filaments embedded in hot am- 
bient plasma has implications for the problem of the 
plasma's thermal conductivity, which must be moder- 
ated by the intracluster magnetic field. If the conductiv- 
ity were high, one would expect the filaments to evapo- 
rate rapidly. Thus, in principle, their very existence limits 
the conductivity in the same way that of cold fronts in 
the ICM does (Markevitch et al. 2000, 2003). Given our 
very limited understanding of the connectedness of the 
intergalactic magnetic field, and the possibility that ther- 
mal conductivity contributes significantly to the global 
dynamics of cooling flows (Narayan & Medvedev 2001, 
Voigt & Fabian 2003 and references therein), observa- 
tional constraints on thermal conductivity are valuable. 

In this paper we address two main questions: 

(i) Why are cold filaments observed only in cooling- 
flow clusters? and why only in the cluster core? 

(ii) To what degree does the existence of these fila- 
ments constrain the thermal conductivity of the ICM? 

The competition between evaporation and condensation 
when a cool body of plasma is immersed in a hotter am- 
bient medium has been studied extensively in the liter- 
ature (Graham & Langer 1973; Cowie & McKee 1977; 
McKee & Cowie 1977; Balbus & McKee 1982; Giuhani 
1984; Draine & Giuhani 1984; Balbus 1986; Bohringer 
& Hartquist 1987; Bohringer & Fabian 1989; McKee & 
Begelman 1990; Ferrara & Shchekinov 1993). Most effort 
has been devoted either to spherical clouds or to pla- 
nar inhomogeneities, although Cowie & Songaila (1977) 
considered spheroidal clouds. Since the morphology of 
observed filaments suggests that cylindrical symmetry 
should provide a useful idealization, we revisit the prob- 
lem for this case, and discuss the implications of our re- 
sults for real systems. 

In any geometry, clouds below some critical size are 
evaporated, while clouds above this size grow by conden- 
sation of the hot medium. We determine the critical size 
that separates these two regimes by finding the steady- 
state solution in which conductivity exactly balances ra- 
diative losses. The typical size of observed filaments is 
likely to be comparable to the critical size: much smaller 
filaments are rapidly evaporated away, while filaments 
above this critical value grow slowly, at a rate determined 
by the cooling time. Significantly larger filaments are ex- 
pected to be very rare: in principle, if they formed they 
could fundamentally change the nature of the system by 
eliminating the cooling hot ambient plasma. 

In Section 2 we derive the critical cylindrical solu- 
tion. In Section 3 we discuss the application of the model 
to cluster filaments. Our conclusions are in Section 4. 



2 RADIATIVELY STABILIZED 
CYLINDRICAL CLOUDS 

2.1 Model and equations 

We consider a long cylindrically symmetric cloud in which 
there is perfect balance between conductive heat flow 
towards the symmetry axis and radiative losses, so the 
plasma flows neither inwards nor outwards (see McKee 
& Cowie 1977). The temperature on the axis is assumed 
to be r ~ 10* K, while far from the origin T asymp- 
totes to a value ~ 10® — 10* K. With these assumptions, 
and neglecting the effects of gravity, the energy equation 
reads 



V-q = n^A(r), 



(1) 



where q is the heat flux, is the electron number den- 
sity, and A(r) is the normalized cooling function. In the 
unsaturated thermal conduction regime^ , the heat flux is 
given by 



-K(r)vr, 



(2) 



with «;(r) thermal conductivity. For a hydro- 
gen plasma, in the absence of magnetic fields, 
Spitzer (1962) derives k(T) = kqT^^^, with 
1.84 X 10"^ (In A) 



ergs ^ cm ^K where 



In A is the Coulomb logarithm, which is only weakly 
dependent on rie and T (in the following we assume 
InA = 30). Confinement of thermal electrons by the in- 
tracluster magnetic field reduces the actual conductivity 
to 



(3) 



where the parameter / < 1 is the heat conduction sup- 
pression factor (e.g., Binney & Cowie 1981, Bohringer & 
Fabian 1989). 

At least in some clusters, observations indicate that 
emission-line nebulae are in pressure equilibrium with 
the surrounding hot plasma (e.g., Heckman et al. 1989). 
So we require that the pressure P = n^T is constant 
throughout the interface, and we use this relation to elim- 
inate He in the right-hand side of equation Q. Thus, 
under the assumption of cylindrical symmetry, eliminat- 
ing q with the help of equations J^J and Q, the energy 
equation Q is reduced to the differential equation 



,A{T) 



(4) 



where R is the cylindrical radius. 
2.2 Boundary conditions 

Equation ||1J applies outside the edge at Rin of the 
T ~ 10* K cloud: we assume that the cloud contains its 
own heat sources, such as photoionizing stars, so on its 
surface the temperature gradient and heat flux vanish. 
We require T to attain some specified value on _R = -Rout . 
Thus we solve equation Q subject to three boundary 
conditions, T(i?in) = Tin = 10* K, dT/dR{Rin) = 0, 
and r(i?out) ~ Tout, the temperature of the ambient 

^ In the Appendix we discuss the possible effects of saturation. 
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Figure 1. Eigenvalues A of equation j^J for a set of values of 
the outer temperature Tout. 



Figure 2. Temperature profile T(R) corresponding to 
the eigenfunction r(g) of equation JSJ, for Tout = 
0.1, 0.2, 0.5, 1, 2, 5, 10, in units of lO'' K. 



hot plasma, which varies with distance from the cluster 
centre. For completeness we will explore the wide range 
10' < Tout < 10« K. 

Introducing the dimensionless quantities (, = 
i?/-Rout, r = r/Tout, and At^^,{t) = A(rrout)/A23, with 



A2 



10" 



'ergs ^ cm^, equation Q reads 



A 



AT„„tM 



where 



A : 



(5) 



(6) 



is an unknown dimensionless parameter that is the 
eigenvalue of equation (|KJ. Expressed in dimensionless 
variables, the boundary conditions are T{^in) = Tin, 
dr/d^{^in) = 0, and r(l) = 1 with = i?in/-Rout and 
Tin = Tin /Tout- Duc to the presence of AT^^t, the eigen- 
functions r(^) and eigenvalues A depend on the outer 
temperature Tout in addition to the dimensionless num- 
bers ^in and Tin. 

The inner and outer radii Rin and iiout cannot be 
exactly identified with observational quantities. However, 
we can assume that Rin corresponds approximately to the 
width of the filament, and i?out to its length I. The latter 
assumption is consistent with the hypothesis of cylindri- 
cal geometry, which does not hold at 7?>/. Observed 
filaments are highly elongated, with ratios length/width 
> 100 (Conselice et al. 2001). Thus, we explore a range 
of values of ^in that satisfy ^in <C 1. 



2.3 Numerical solutions 

Equation is non-linear and would require numerical 
integration even if ATo^t (''') were approximated by a sim- 
ple analytical function. Thus, use of the tabulated cooling 



function of Sutherland & Dopita (1993) for metallicity 
[Fe/H] = —0.5 introduces no additional computational 
complexity. We use a shooting method, with Newton it- 
eration, to calculate the solutions. 

Numerical integration of equation ^ shows that, for 
any fixed Tout in the range considered, the value of A is 
virtually independent of ^in, provided that ^in is suffi- 
ciently small (i.e., ^in^lO"^). Thus, in what follows we 
can simplify our treatment by assuming that the eigen- 
value A depends only on Tout. 

Although the details of the inner temperature pro- 
file r(^) depend significantly on ^in, the behaviour of the 
profile at ^ > 10"'^ is virtually independent of ^in pro- 
vided Cin^lO"^. Henceforth we set Cm ^ 10"^ so that 
both A and the observationally significant part o/r(C) are 
independent on the exact radial range considered. 

Fig. shows A as a function of Tout, while Fig. |21 
shows a number of the corresponding eigenfunctions 
T{R). For values of Tout in the relevant range, A spans a 
relattvely small range (0.08 < A < 0.37), which is reduced 
to 0.25 < A < 0.37 for Tout > 4 x lO'' K. 

From Fig. |5|is apparent that over a significant frac- 
tion of the radial range the temperature profiles are 
roughly power laws, with similar indices. This suggests 
that an approximate analytical solution could be ob- 
tained by exploring the asymptotic behaviour of the 
eigenfunctions of equation for ^ ^ 1 (r ^ 1). When 
we assume r = C", and approximate the cooling function 
near Tout by ATo„t(r) ~ [A(rout)/A23]r°', equation @ re- 
quires?! = 4/(ll-2a), and A = 5n^/[A(Tout)/A23]. Since 
— l/2<Q!<l/2 for the relevant temperatures, the outer 
slope of the eigenfunction is expected to lie in the small 
range l/3<n<2/5 (cf. Fig.|5J. The run of the eigenvalue 
as a function of Tout (Fig. reflects the proportionality 
to l/A(Tout), modulated by n^, which accounts for the 
dependence on the local slope of the cooling function a. 
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Figure 3. Left: Critical length of a radiatively stabilized fila- 
ment as a function of the temperature of the ambient medium, 
at fixed pressure P = lO^Kcm"^, for a few values of the 
suppression factor /. The dot-dashed line represents the self- 
similar relation between virial temperature and virial radius. 
Right: Critical length as a function of pressure, at fixed tem- 
perature T = 3 X 10^ K, for the same set of values of / as in 
left panel. 



shows Zcrit(7out) from equation ||7|l, for pressure P — 
10** K cm , which is a typical value in the core of a 
cooling-flow galaxy cluster [e.g., for the Perseus clus- 
ter, Conselice et al. (2001), Schmidt, Fabian & Sanders 
(2002); see also Heckman et al. (1989)]. A few values of 
the suppression factor (/ = 1, 0.1, 0.01, 0.001) are consid- 
ered. For given /, the relation Zcrit(rout) divides the the 
diagram in two regions: filaments longer than ?crit are 
condensing, while smaller ones are evaporating. It is ap- 
parent from the steep slope of the curves in the diagram 
that the process of evaporation or condensation of fila- 
ments is very sensitive to the temperature of the ambient 
medium. For a plausible value of the suppression factor, 
/ ~ 0.01, at fixed pressure ^ lO^Kcm""^, any filament 
longer than '--^ 0.1 kpc is condensing in a cluster with am- 
bient temperature Tout = 10^ K, while filaments as long 
as ~ 10 kpc are evaporating if Tout = 6 x 10^ K. 

Although /crit depends on pressure more weakly than 
on temperature, pressure also plays an important role in 
determining the fate of filaments. This is apparent in the 
right diagram of Fig. |21 where we fix Tout = 3 x 10^ 
and plot Zcrit as a function of P, for the set of val- 
ues of / considered above. The diagram shows that, 
at fixed temperature, in low-pressure clusters virtually 
all filaments are evaporating (e.g., Zcrit ~ 100 kpc for 
P ~ few X 10* Kcm~^ and / — 0.01), but condensation is 
effective in the highest-pressure ambient gas (Zcrit ~ 1 kpc 
for P>10'^Kcm-=' and / = 0.01). 



3 APPLICATIONS 

3.1 Critical length for filaments 

We now discuss the physical implications of the results 
presented in the previous section. Equation Q relates 
the numerical value of A to the dimensional quantities 
associated with a filament and its surroundings. If, as 
explained above, we identify i?out with the length of a fil- 
ament, then this equation yields for a given environment 
the critical length Zcrit of filaments, such that smaller fila- 
ments evaporate while larger ones grow. Numerically we 
have 

^crit = Pout 

-'■V"=A'"(^)""(j;j^)->,pc,7, 

This equation indicates that Zcrit has a strong dependence 
on the outer temperature {lait oc T^^/*, neglecting the 
weak dependence of A on Tout; see Fig.0, and also is in- 
versely proportional to the pressure P. Even though the 
dependence of l^it on the thermal conduction suppression 
factor / is not very strong {lait oc f^^"^), f can signifi- 
cantly affect Zcrit because the value of / is very uncertain: 
/ ~ 1 seems to be excluded, but different authors gives 
values differing by two orders of magnitude, from a few 
10"^ to - 10"^ - see Narayan & Medvedev (2001) for a 
discussion. We consider values in the range 10"'' / 1 
and discuss the possibility of constraining / from obser- 
vations of filaments in Section 3.4. 

Consider the implications of the dependence of ?crit 
on Tout at fixed pressure. The left panel of Fig. |21 



3.2 Time-scales 

What are the characteristic time-scales for evaporation 
and condensation of filaments? A detailed answer to this 
question would require the solution of the full hydrody- 
namic equations for non-vanishing mass flow. However, 
for the purpose of our investigation we are interested in 
an order of magnitude estimate, which can be obtained 
by considering that the time variation of temperature is 
described by a differential equation of the form 

dT T T 
^ 7 7 — ' 

UjL iti t'COol 

where t^, and tcooi are the thermal conduction and cool- 
ing time-scales, respectively. The relation between these 
time-scales is given by 

fcooi |V-q| ~ k{T)T ' ^ 

because |V ■ q[ ~ k{T)T/1^, where I is the temperature 
variation length scale (Begelman & McKee 1990). In the 
considered case, I is the radial extension of the interface, 
which is of the order of the fllament length (see Section 2). 
Considering that the equilibrium solution (/ — la-it) cor- 
responds to ~ tcooi, from equation JUJ we get 




and equation ^ can be approximated as 
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Figure 4. Dotted curves: density— temperature relation for 
plasma with fixed cooling time (tcool = 0.1, 0.3, 1, 3 in units of 
10^" yr). Solid curves: density as function of temperature for 
'crit = 1, 10, 100 kpc and / = 0.01 (equation^. 

Consequently, filaments shorter than ~ |icrit evaporate 
in a fraction (//Zcrit)'^ of the cooling time, which for ex- 
ample is < 300 Myr within ~ 4 kpc of the centre of the 
Hydra cluster, (David et al. 2000) . Filaments smaller than 
'crit by a factor > 8 have evaporation times < 5 Myr that 
are smaller than their likely dynamical times, Ijv^ where 
V ~ lOOkms"^ is the typical gas velocity estimated from 
line widths (Sabra et al. 2000 and references therein). On 
the other hand, equation \\\\ in the limit I ^ icrit indi- 
cates that, as expected, very long filaments condense on 
time-scales of the order of the cooling time of the ambient 
medium. 

These considerations strongly support the hypothe- 
sis that only filaments with lengths comparable to the crit- 
ical length or longer are likely to survive for long enough 
to be observed. 

3.3 The fate of filaments in cooling-flow and 
non— cooling-flow clusters 

The results presented above have interesting implica- 
tions, when compared with the properties of observed 
galaxy clusters. 

(i) A first simplistic implication is that, due to the 
strong dependence of lait on temperature, in massive 
(high-temperature) clusters most cold gas structures are 
evaporating, while even relatively small clouds can sur- 
vive (and condense) in small (low-temperature) clusters. 
On the other hand one needs to bear in mind that clus- 
ters with higher temperatures have larger virial radii Rv, 
and it is natural to expect the typical size of clouds to 
increase proportional to Ry. The left panel of Fig. |3| 
shows the relation between Rv and the virial temper- 
ature (i?v oc Ty^) predicted for galaxy clusters un- 



der the assumption of self-similarity (the slope of this 
scaling relation is consistent with that derived from X- 
ray observations of galaxy clusters; see, e.g., Ettori, De 
Grandi & Molendi 2002, and references therein). Thus, 
as is apparent from the diagram, the relation icrit(T'out) 
is much steeper than the relation between cluster size 
and temperature. Hence, in the highest-temperature clus- 
ters any cold filament is evaporating, while in the lowest- 
temperature clusters, any resolvable filament is condens- 
ing. 

(ii) The scenario presented in the previous point would 
hold under the assumption that the physical properties 
(gas density and pressure) of clusters of given tempera- 
ture would not differ significantly from cluster to cluster. 
In fact, we know this is not the case: for fixed gas tem- 
perature, non-cooling-fiow clusters (i.e., ones with long 
cooling times) have, by definition, lower central density 
(and pressure) than cooling-fiow clusters. Thus, as shown 
by the right panel of Fig. |3] since lait oc l/P, a filament 
long enough to be condensing in the core of a cooling- 
fiow cluster with central temperature T might evaporate 
if it were at the centre of a non-cooling-fiow cluster with 
the same temperature (but lower density). This argument 
suggests that filaments of cold gas are not observed in 
non-cooling-flow clusters because in these clusters they 
are evaporated very efficiently, even at the cluster cen- 
tre. On the other hand, whatever their formation mech- 
anism, filaments are more likely to survive and grow in 
the denser, higher pressure cores of cooling-flow clusters. 

(iii) In a given cooling-flow cluster, T increases as one 
moves away from the centre, while P falls. Consequently, 
equation lO predicts that icrit increases rapidly with dis- 
tance from the cluster centre. This result explains why, 
even in the infall or merging scenario for the formation 
of the filaments, these are observed only in the cooling 
cores of clusters, and not in the outer regions. For ex- 
ample, in the Perseus cluster, the gas temperature at a 
distance ^ 100 kpc from the center is roughly twice the 
central temperature. Also, in the same radial range the 
pressure decreases by a factor of ~ 2.5 (e.g., Schimdt et 
al. 2002; Churazov et al. 2003). Thus, equation ^ pre- 
dicts that icrit grows by over a factor of 10 between the 
centre and 100 kpc. Interestingly, even in the core, obser- 
vations indicate that the mean Hq emission is a strongly 
decreasing function of the cluster radius (e.g., Conselice 
et al. 2001). This result may arise because as one moves 
out, Hq emission from smaller filaments contributes less 
and less to the overall signal. 

We can summarize the points above by considering 
the relation between lait and the cooling time of ambient 
gas. For simplicity, we use the approximate isobaric cool- 
ing time fcooi — 10*n~^r^''^ yr (where is in cm^"* and 
T in K; e.g., Sarazin 1986), considering the temperature 
range 10^ — 10* K. The dotted curves in Fig. 01 represent 
the loci of constant cooling time in density-temperature 
space, while the solid curves show the loci of constant 
'crit, for suppression factor / — 0.01. We see that if the 
cooling time is long (as in the cores of non cooling-flow 
clusters, and in the outer regions of any cluster), lait is 
large, so evaporation is very likely. On the other and, 
small Zcrit corresponds to a short cooling time. Hence it 
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Ti„=10-* K 
P = 106 K cm-3 - 
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Figure 5. Differential emission measure EM(T) = 
n^dV/dlogT as a function of T, corresponding to the critical 
stationary solution for outer temperatures Tout = 1,2,5,10, 
in units of 10^ K. / measures the reduction of heat conduction 
below Spitzer's value. /3 is the index of the best-fitting power 
law. 



is precisely the physical conditions found in the cores of 
cooling-flow clusters that favour condensation. 



3.4 Suppression of thermal conductivity in the 
intracluster medium 

As outlined in the Introduction, the existence of cold fil- 
aments could be used to constrain the value of the ther- 
mal conductivity in the ICM. Indeed, we could derive 
the value of the suppression factor / from equation (|7|l 
if we knew all the other quantities involved. Unfortu- 
nately, while Tout and P are well constrained by obser- 
vations, it is quite difficult to put strong observational 
constraints on lait- A detailed morphological study of 
the filaments, available only for the Perseus cluster (Con- 
selice et al. 2001), reveals complex structures on different 
scales. However, we note that even a quite conservative 
estimate of lait can lead to interesting constraints on /. 
For instance, the data suggest that in the Perseus clus- 
ter filaments ~ lOkpc long are present at least up to 
radii ~ 40 — 50kpc. This means Zcrit ;$ lOkpc, where the 
ICM temperature and pressure are ~ 4.5 x 10^ K and 
~ lO^Kcm"^, respectively (Schmidt et al. 2002). Ac- 
cording to equation J7|, for these physical properties of 
the ambient medium, this condition is translated into the 
upper limit f < 0.04 for the thermal conduction suppres- 
sion factor. This result is compatible with that found by 
Markevitch et al. (2003) on the basis of observations of 
cold fronts in A754. We note that a substantially higher 
efficiency would be required in order for the thermal con- 
ductivity to significantly affect cooling flows in the cen- 
tres of clusters (e.g., Narayan & Medvedev 2001). 

In the discussion above we considered the thermal 



conduction suppression factor as a global parameter char- 
acterizing the ICM in clusters. In fact, the reduction of 
thermal conductivity could vary from cluster to cluster, 
or within a given cluster, depending on the detailed prop- 
erties of the magnetic field. However, there are indications 
that the effect of the reconnection of the magnetic field 
lines between the filaments and the ICM would be to 
enhance the efficiency of heat conduction in the centres 
of cooling-fiow clusters (Soker, Blanton & Sarazin 2004). 
Thus, at least according to this scenario, it seems reason- 
able to extend the derived upper limit on / to the ICM 
in general. 



3.5 X-ray emission from conduction fronts of 
filaments 

Here we explore the possibility that the conduction fronts 
described above can be detected in soft X-rays. This is 
particularly interesting since Fabian et al. (2003) find soft 
X-ray emission associated with the brightest Hq emis- 
sion filaments around NGC1275, suggesting that it is pro- 
duced in the conductive evaporation or condensation of 
the filaments. We quantify the emission properties of the 
conduction fronts for the considered cylindrically sym- 
metric clouds, computing the differential emission mea- 
sure 



EM(T) 



, dV 
'dlogT' 



(12) 



where is the electron number density and dV is the 
volume element. Under the hypothesis of constant pres- 
sure P = UeT, the differential emission measure can be 
related directly to the temperature profile T{R): using 
dV = 2nlRdR, with I length of the filament, we get 



EM{T) = 2ttIP- 



R dR 
dlogT' 



(13) 



As discussed in Section 2, the temperature profile, for 
given Tout, can be considered independent of ^in at large 
enough radii (^ > lO"'^); correspondingly we find that the 
differential emission measure at T>0.1Tout is not sensi- 
tive to the value of ^in. Fig.|5]plots the emission measure 
as a function of T for a subset of values of the outer 
temperature (Tout = 1, 2, 5, 10 in units of 10^ K). For the 
normalization of EM{T) we assume P — lO^Kcm"^, 
and l = Icrit = Rout, in accordance with the discussion 
in Section 2. Due to the dependence of lait on the ther- 
mal conduction suppression factor / (equation ITJ, the 
normalization of EM{T) is proportional to f^^^ ■ Thus 
it is difficult to make predictions on the absolute value 
of the emission measure, also because of the strong de- 
pendence on pressure. On the other hand, the shape of 
EM{T) is not affected by these parameters, being deter- 
mined by the temperature profile, which depends only 
on the temperature of the ambient medium. As apparent 
from Fig. |^ in the considered temperature range the dif- 
ferential emission measure can be roughly approximated 
with a power-law EM{T) oc T''. We find that EM (J) is 
steeper for lower Tout, with best-fitting slope /3 ranging 
from P ~ 3.4 for Tout = 10* K to /3 ~ 4.6 for Tout = 10^ K. 
Accurate measurements of the soft X-ray spectra around 
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filaments could be interestingly compared with the slopes 
of these emission measure profiles. 



4 CONCLUSIONS 

Embedded within cooling flows, one frequently detects 
filaments of cold gas. Thermal conduction will cause heat 
to flow into these filaments from the hot ambient gas, 
with the result that filaments below a critical size will 
be evaporated. On the other hand radiative cooling is 
very efficient in cold, dense gas, so larger filaments can 
survive, and indeed nucleate the cooling of the ambient 
medium. We have determined the critical filamentary size 
that divides the two regimes. 

The absolute value of this length scale, lait, is un- 
certain because we do not know how effectively the mag- 
netic field suppresses thermal conduction. By contrast, 
the way in which Zcrit scales with the temperature and 
pressure of the ambient medium is well determined and 
leads to suggestive results. 

The critical length scale rises sharply with ambient 
temperature: lait (x T^^^'^, and decreases as with 
ambient pressure. Consequently, in the Perseus cluster 
Icvit increases by over an order of magnitude between the 
cluster centre and a radius of 100 kpc, where it is probably 
larger than 10 kpc. Within a few kiloparsecs of the cluster 
centre, filaments only a kiloparsec long may either grow, 
or at least take a long time to evaporate because their 
lengths are just below the critical value, while such small 
filaments are rapidly evaporated in the main body of the 
cluster. 

Where do filaments come from in the first place? 
They cannot emerge spontaneously from the ambient 
medium because this is thermally stable,'^ so they must be 
injected. Simulations of structure formation predict that 
significant amounts of cold gas fall into clusters (Katz et 
al. 2003) , and the observed filaments could represent just 
the fraction of this infall that reaches the centre. There is 
also some indication that extremely close to the cluster 
centre a small amount of gas cools to filamentary tem- 
peratures, and some of this gas could be ejected by the 
AGN to the radii where filaments are observed. 

Whether infall or AGN are responsible for injecting 
filamentary gas, there is no reason to expect the lengths 
of the injected filaments to have a characteristic value. 
It is natural to think in terms of a spectrum of filamen- 
tary lengths I, in which the number of filaments decreases 
with increasing filament length, probably as a power law. 
Observations of filaments would then be dominated by 
filaments with I ~ icrit because shorter filaments would 
be rapidly evaporated, and larger ones would be rare. The 
local number density of filaments would be a decreasing 

^ Balbus (1991) shows that thermal instability is possible in 
the presence of even a weak magnetic field. However, a weak 
field can only control low-amplitude perturbations, and buoy- 
ant stabilization will reassert itself once the perturbation am- 
plitude is greater than Pmag/^gas ~ 10~^ in clusters. It is also 
the case that for technical reasons Balbus had to insist on a 
highly artificial background field for which field-line closure is 
problematic. 



function of lait, explaining why observed filaments are 
concentrated near the centres of cooling-flow clusters. 

Since filament growth depends on the efficiency of 
radiative cooling, Zcrit tends to be small when the cooling 
time of the ambient medium is short. Hence clusters with 
short central cooling times, have smaller values of Zcrit and 
more filaments than clusters with similar temperatures 
but longer cooling times. Consequently, it is natural that 
filaments are only detected in cooling-flow clusters. 

For the case of a filament of critical length, we have 
calculated the temperature profile throughout the inter- 
face between gas at ~ 10* K that emits optical emission 
lines, and the ambient X-ray emitting plasma. This pro- 
file yields the run of differential emission measure with 
temperature, which can be tested observationally. 

As cosmic clustering proceeds, the temperature of 
gas that is trapped in gravitational potential wells rises 
as longer and longer length-scale perturbations collapse. 
This temperature rise causes lait to rise much faster than 
the virial radius of the potential well. Hence, while at 
early times filaments that are much smaller than the sys- 
tem can grow, growth is later restricted to filaments that 
are comparable in size to the entire system. Consequently, 
the ability of infalling cold gas to nucleate star formation 
declines as clustering proceeds. 

Semi-analytic models of galaxy formation that can 
account for the observed slope of the galaxy luminosity 
function at the faint end need efficient feedback from star 
formation, and over-produce luminous galaxies (Benson 
et al. 2003). Binney (2003) has argued that galaxy for- 
mation is dominated by infall of cold gas, and that it is 
the efficient evaporation of filaments in massive systems 
that suppresses the formation of massive galaxies. Our 
results make a step towards quantifying this picture. 
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g.at=0.4(^^) nefesT, (Al) 

where ks is the Boltzmann constant, and me is the elec- 
tron mass. Cowie & McKee introduced a discriminant of 
whether the heat flux in a plasma is saturated: this is the 
saturation parameter a = qd/qsnA, where ^ci = ^^{T)\^T\ 
is the modulus of the classical heat flux (equation|5J. The 
critical value that separates the two regimes is o" ~ 1: for 
(T <g 1 the heat flux is unsaturated, while for a ^ 1, 
equation lAll l applies. 

For given Tout and corresponding critical tempera- 
ture profile T{R), the value of the saturation parame- 
ter as a function radius in the interface can be obtained 
combining equations (|2J, ©, and (lAlll . and using the 
relations Tie = P/T and flout = 'crit (equation |7||. Nu- 
merically we get 

aiR) ^ 8 X 10-//^A-/^ (^)''%^|, (A2) 

where ^ = R/Rout, t = t(£^) is the normalized temper- 
ature profile, and A = A(rout) the corresponding eigen- 
value (see Section 2). We note that o{R) is independent 
of the pressure and scales as the square root of the sup- 
pression factor /. Applying equation IIA2t to conduc- 
tion fronts with ambient gas temperature in the con- 
sidered range 10' < Tout < 10** K, we find in all cases 
""(-R) ~ lO"'^/^''^- Thus, the condition for unsaturated 
heat conduction (a <C 1) is satisfied throughout the inter- 
face, even for the largest values of the suppression factor 
/~1. 

This paper has been typeset from a I^^X/ KT^^X file pre- 
pared by the author. 



APPENDIX A: THERMAL CONDUCTION 
SATURATION PARAMETER 

We have hitherto assumed that the thermal conduction 
is not saturated, so the heat flux is given by the classi- 
cal equation with the thermal conductivity given by 
equation Q. The latter formula is derived under the as- 
sumption that le/lr -C 1, where 1^ is the electron mean 
free path and It ~ r/|Vr| is the temperature scale 
height. When this condition is violated, equation Q over- 
estimates the heat fiux, which is limited by the thermal 
speed of electrons, and the modulus of q is given by the 
saturated heat flux (Cowie & McKee 1977) 
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